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ABSTRACT 

Galactic disks in triaxial dark matter halos become deformed by the elliptical potential in the 
plane of the disk in such a way as to counteract the halo ellipticity. We develop a technique to 
calculate the equilibriu m configura tion of such a disk in the combined disk-halo potential, which is 
based on the method of iJod (|200Gf ) but accounts for the radial variation in both the halo potential 
and the disk ellipticity. This crucial ingredient results in qualitatively different behavior of the disk: 
the disk circularizes the potential at small radii, even for a reasonably low disk mass. This effect has 
important implications for proposals to reconcile cuspy halo density profiles with low surface brightness 
galaxy rotation curves using halo triaxiality. The disk ellipticitics in our models are consistent with 
observational estimates based on two-dimensional velocity fields and isophotal axis ratios. 
Subject headings: galaxies: kinematics and dynamics — galaxies: halos — galaxies: spiral — galaxies: 
structure — methods: numerical — dark matter 
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r _~ 1 1. INTRODUCTION 

< 

Galaxies are thought to be surrounded by large dark matter halos. These halos are much more massive than the 
visible components of galaxies, and dominate much of the dynamics. Although dark matter halos are often assumed to 
be spherical for simplicity, the halos that form in cosmological simulations are quite flattened, with typical intermediate 
. axis ratios of b/a ~ 0.8 and minor axis ratios of c/a ~ 0.6 . with some systematic variation depending on the mass of 
the halo and the radius at which the shape is measured ([Warren et al J [19921 ; Uing fc Sutol 120021 iBailin fc Steinmet3 
l-L; [20051: lAlkood et al.ll2f)0ll . This non-sphericity is a testable prediction of cosmological models. 

Simulations of disk galaxy formation within dark matter halos fi n d that the presence of the disk modifies 
the shape of the halo, red ucing the halo triaxiality ([Dubinskil Il994t iKazantzidis et all [2004: Baili n et al.l [2005; 
iBerentzen fc Sh losman 2006). However, as long as the final shape of the halo retains some ellipticity in the plane of the 
disk, the dynamics and shape o f the disk will be affected by the deviations from axisymmetry (e.g. iGerhard fc Vietrl 
■ fl9M ISchoenmakers et allll997h . 

Observations indicate that many disks do indeed have small but non-zero ellipti cities. Evidence for elliptical disks 
. . comes from harmonic deco mposition of galaxy photometry plix fc Zarit skv 1995), harmonic decomposition of two- 
^sO ' dimensional velocity fields (ISchoenmakers et ailll997t ISimon et al.ll2005l ). and statistical analysis of the distribution 
CD of projected shapes (|Rvdenll2006[ ). These results qualitatively confirm that galactic dark matter halos are elliptical; 
precision m easurements could provide direct constraints on the shapes of the halos. 
Recently, lHavashi fc Navarrol ([2006L hereafter HN06) proposed that elliptical orbits within the disk produced by the 
\ triaxiality of the halo could reconcile cuspy density profiles that form in cosmological simulations ([Navarro etaII[l9M 
hereafter NFW) with observed rotation curves of low surfac e brightness (LSB) galaxies, which often appear to require 
halos with constant density cores (e.g. Jde Blok et all llOOl). This analysis did not take into account the self-gravity 
J^j , of the disk. In galaxies with massive disks, the gravity of the disk contributes to the net potential and the dynamics 
. - - i of the disk are determined by a combination of the halo and the disk itself. In order to draw conclusions about the 
shape of the halo from the measured dynamics of the disk, we must determine self-consistently both how the disk is 
perturbed by the potential and how the perturbed disk contributes to the potenti al. 

An elegant method to carry out these calculations was proposed by Pod (|2000D (hereafter |J2K| ; see also Uodll997l 
1999). By assuming a logarithmic halo p otent ial with a small constant elliptical perturbation and an exponential disk 
with a small constant elliptical response, IJ2KI solved for the self-consistent response. She demonstrated that the disk 
response dilutes the ellipticity of the potential mos t stro ngly at 1.42 disk scale lengths. 

There are a number of simplifying assumptions in lJ2Kl that require examination. The mos t impor tant assumption is 
that both the halo perturbation and the disk response are constant with radius. In contrast. IHN06I demonstrated that 
a radially-varying perturbation is required to reconcile LSB long-slit rotation curves with cuspy halo profiles. Indeed, 
cosmological simulations predict a radially-varying perturbation in the halo potentia l; even if halo s had isodensity 
surfaces of constant ellipticity, the shape of the potential would vary with radius, and lHavashi et al. (2007, hereafter 
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HNS), who directly measured the shapes of isopotential surfaces of cosmological halos, found even stronger variation 
with radius. The response of the disk is also not expected to be uniform in a radially-varying potential. 

In this paper, we generalize the method of IJ2K1 to the more realistic case of radially varying halo perturbations and 
radially varying disk responses. In § [2] we detail the method for determining the disk shape and dynamics. In § [3] 
we use this method to determine the shapes and dynamics of disks in sample triaxial halos and demonstrate how the 
results depend on the properties of the disks and halos. § H discusses our results in the context of observations that 
directly probe disk ellipticity, and in §[5] we present our conclusions. 

2. METHOD 
2.1. Outline 

Our method, which is based closely on[j2K, is as follows: 

1. Calculate the axisymmetric component of the potential and the elliptical perturbation in the potential induced 
by the triaxial halo (§ 



2. Calculate the closed orbits and corresponding disk ellipticity for a given net perturbation to the potential (§§ 12.31 
aii(irO). 



3. Calculate the elliptical perturbation in the potential induced by a given disk ellipticity (§ \2.5\i . 

4. Solve for the form of the net potential perturbation that satisfies all of the above constraints (§ 12.6 



Throughout this procedure, the halo is kept fixed, i.e. it does not respond to the presence of the disk. Therefore, the 
potential that should be used in the calculation is the real shape of the halo in which the disk lies, which is less triaxial 
than the shape the halo wou ld have in the absence of baryonic processes (jKazantzidis et all [2004; Ba ilin et al.l 120051 : 
iBerentzen fc Shlosmanll2006[ K 

The main differences between this work and IJ2KI are: 

1. Where QUI assumes a logarithmic potential for the halo, we evaluate the radial form of the potential directly 
from a density distribution motivated by cosmological simulations. 

2. Where EH assumes a constant perturbation to the halo potential, we allow the perturbation to vary with radius 
and either evaluate it directly from a triaxial density distribution motivated by cosmological simulations or use 
parametrizations developed from measurements of halos in cosmological simulations. 

3. Where EH assumes that the disk responds with a constant ellipticity, we allow the disk ellipticity to vary with 
radius. 

Whenever we carry out numerical calculations in this paper, all radially- varying functions are tabulated on a radial 
grid sampled at 100 radii Ri spaced logarithmically between 0.1 and 100 kpc in order to finely sample the inner 
region of the disk where the quantities vary most rapidly; 50 grid points lie at R < 1 kpc. The functions are linearly 
interpolated between grid points when their values are required at arbitrary radii. 

2.2. Axisymmetric potential and halo perturbation 

For the perturbative approach, we assume that within the plane of the disk (which we take to be z = for simplicity), 
the total potential can be written as 

0)=$oOR)+$pcrtOM) 

= $ (R)(l + f peit (R)cosm6), (1) 

where R, 9, and z are the cylindrical coordinates. For the purposes of this paper, we will assume an elliptical 
perturbation (i.e., m = 2) from now on (see Appendix [Al for a discussion of the m = 4 mode). We assume that 
/pert(-R) is small and varies slowly with R. Note that for / pe rt > 0, isopotentials are elongated along the x axis and 
closed orbits in the disk are elongated along the y axis. 
Both the disk and halo contribute to both the axisymmetric and m = 2 components of the potential: 

$ (i?) = Q%*>(R) + < isk (i?), (2) 

/port(-R) = fpelt(R) + /pcrt(-R)- ( 3 ) 

To first order, the m = 2 perturbation in the potential induces an m — 2 perturbation in the surface density 
distribution of the otherwise exponential disk (see Appendix [Al for a justification of our decision to neglect the higher- 
order terms): 

■K 1 - 5 ^-)] (4) 



T,{R,6) = S exp 



We assume edisk(-R), the ellipticity of the isodensity ellipse, is small and varies slowly with R. The axisymmetric 
component of the disk potential is given by 

$f k (i?) = -nGZoRMy^iy) - h(y)K (y)} (5) 
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(y = R/2Rd] see lFreemanlll970t iBinnev fc Tremainel[T987l eq. 2-168), where / and K are modified Bessel functions. 

Given the density distribution of the dark matter halo, both axisymmetric and m — 2 components of the halo 
potential can be evaluated. If the isodensity surfaces of the halo are self-similar ellipsoids, then the halo density can 
be written as: 

p(x,y,z) = p(s), (6) 

where 

'-(=)"+ «)'+(!)'■ m 

For example, we can use an NFW form for the density: 

P( s ) = , P ° ,2 ( 8 ) 
r- + 



(N avarro et al.lll996l : iJing fe Sutoll2002h . We make no assumptions about the relative magnitudes of a, b, and c; they 
are simply the relative axis ratios along the x, y, and z axes respectively. Therefore, the disk, which lies in the xy 
plane, can be oriented in any of the principal planes of the halo. 
We calculate the halo potential along the x axis, $^ alo (i?) = $ halo (i?, 0, 0), and along the y axis, $J) alo (i?) = 

$ halo (0, R, 0), by numerically integrating eq. (2-99) of iBinney fe Tremaind (|1987l ) using p(s) given in equation ([8]). 
This allows us to calculate $fj al ° and f^rt as 

^°(R) = l[$ x (R) + $ v (R)}, (9) 

and 

f ^ {R) = 2^{R) [q>x{R) ~ ® y{R)] ■ (10) 
2.3. Closed orbits 

In a triaxial potential, dissipative gas settles on stable closed loop orbits when such orbits exist (|E1-Zantll2001[ ). This 
is the case throughout the potential of a centrally-concentrated mass profile such as the NFW profile. Therefore, the 
structure of a galactic disk, which consists of gas clouds and stars formed within those g as cl ouds, is determined by 
the form of the clos ed orb its. These have been examined in detail bv lSchoenmakers et al.l (|1997l ). and simplified into a 
convenient form bv lHN06l . These previous derivations have assumed that the perturbation to the potential is constant 
over the radial excursion of an orbit, which is not the case for the radially- varying perturbations that we wish to study. 
We have therefore rederived the equations for closed orbits within a radially-varying perturbation from the equations 
of motion. The orbits follow 

R = R c (l-hE^l mB 29 ) (11) 



with velocities 



e = 8 + ai2 + a32 / P e rt sin 20 o , (12) 
2m 

Vr = V c /pert ai2 sin 29 (13) 

Ve = v J 1 + fpM cos26o \ (w) 



where Rq and 0q define the guiding center of the orbit, 9q = Hot, and the following functions of $o(R) are evaluated 
at R n : 

V C (R) = RQ (R) (16) 



V esc (R) = v/2|$^ (17) 
^ = ^- (m2 " 3) (18) 



dlm(R) = 



2 



9m(R) 



(, . 1 R d/pert \ 

V? \ + 2 /pe rt dR J 



(19) 



a 3m (R) = a lm (R) + . (20) 
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The coefficients a\ m and a^m quantify the degree to which the radius R and angular velocity Vg respectively, wh ich are 
constant for a circular orbit, vary for a unit perturbation to the potential. Our expressions differ from those in HN06 
for the following reasons: (1) we have taken the radial variation of the perturbation into account in our deri vation of 
equation (j 19[) ; (2) we have generalized the expression for g m (R) to be valid for all m, while the expression in HN06 is 
specific to m = 2; (3) we have removed the factor of / pcr t from the definitions of a\ m and a^m for convenience later 
(note that these quantities still depend implicitly on / per t through its derivat ive); a nd (4) we have generalized the 
expression for a^ m to be valid for all potential profiles, while the expression in IHN061 is only valid when g 2 (R) = — 3 5 
which is not the case in the inner regions of an NFW potential. 

Given the tabulated values of <I>o, these functions can be evaluated at the grid points Ri. Because $o has been 
calculated from analytic functions, the tabulated values are relatively free of noise and even the numerical second 
derivative does not contain large fluctuations. Note that the closed orbits are elongated perpendicular to the isopotential 
contours. 

2.4. Disk ellipticity 
The disk must satisfy the continuity equation. In cylindrical coordinates: 

^[R-Z(R,6)V R (R,e)} + -^[Il(R,8)V (R,d)} = O. (21) 

We substitute S(i?, 9) from g]), Vr(R, 9) from ([j"3"]l. and V e (R, 9) from HU). To first order in the small quantities / p0 rt, 
£disk, and their derivatives: 



R 

^-£disk(i?) = /pcrt(-R) 
rid 



. R R dV r , 

a ^ R ^ l -R- d + V c ^- a ^ R) 



(22) 



The neglected second-order terms induce small m = 4 perturbations in the disk; see Appendix fAl for details. 6 Equa- 
tion (|22[) provides us with a relationship between the radial profile of the potential (embodied in a\ 2l 032 and V c ), the 
strength of the perturbation in the net potential (/ pe rt), and the ellipticity of the disk (edisk)- 

2.5. Disk perturbation potential 
The ellipticity of the disk generates an m = 2 perturbation to the disk potential. We calculate this as follows: 

00 />oc 

$ dWk _ $ disk _ $ disk = -G exp(im0) / J m (kR) exp(-fc \z\)dk 

m— — oo 

\ J m (kR')R'dR' [E(R', 9') - £ exp(-i?'/i? d )] exp(-im0') AO 1 (23) 
Jo Jo 

(jBinnev fc Tremainelll987l . eq. 2P-8). We restrict ourselves to the plane = 0, corresponding to an infinitely thin 
disk. 

For small perturbations, the perturbed surface density is 



£« 9') - So e- R '' R * « S e' R ' IRd ed ^ R>) cos 29'. 

Rd 2 



(24) 



Substituting ([24| into (|23|) . we note that the integral over d9' vanishes except when m = ±2. Since 
J 27r cos 26»' e ±ae ' d9' = tt, e lx + e~ lx = 2cosa;, and J 2 (x) = J- 2 (x), we find 

$^ t (R,d) = -7rGE cos2d j J 2 (kR)dk j J 2 {kR')R' exp ( -— J — e disk (i?') dR 1 . (25) 

We can express (R' / Rd) edisk(-R') in the final integral in terms of / pcr t, 012, 032, and V c using 

Z" 00 f°° ( fi" \ 

$f c f t {R,9) = -wGE cos2dJ J 2 (kR)dkJ J 2 (kR') R' exp (-—J 



/pert(-R') 



, R' , R' dV c {R') ■ 
ai2{R) ^-R- d + V c {R) dR )- fl 32(i?) 



dR'. (26) 



Because we do not know a priori the net perturbation / pert (i?), we cannot immediately evaluate these integrals. 
However, if we can find a function / P c° t xy whose form is similar to / por t, i-e. if /pert/Xpert^ * s a slowly- varying function 
of R, then we can approximate the potential as: 



$$&{R,0) « - ^ R K gZ cos29 f J 2 (kR)dk f J 2 {kR')R' exp (-?f) 

/port I 7 *-) JO JO \ U dJ 

dR'. (27) 



/ pert 

rpr 
1 pcrt 

We have also omitted the term in equation {22) that is proportional to d(/ pC rt 0*12) /dR; however, its effect is negligible 
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Note that ai2(i?') and a^iR') depend implicitly on fpc°t Y (R') through its derivative. A first approximation can be 
obtained by setting fp™* y — /pert > whose values have been tabulated from (JTUJ). Because the integral over R' is 
independent of R and the integral over k is independent of R' , the integrals can be evaluated independently on fine 
grids of k and R, respectively. Using this technique, we calculate 



f°° f 00 ( R! \ 

I J 2 (kR)dkj MkR')R'ex P i- w ) 



dR' (28) 



CT) [*»(*) (i - 1 + ^^f) - 

at each grid point Ri. The perturbation potential due to the disk is then given by 

$^(R, 9) = -/pert(iJ) V(R) cos 29. (29) 

Expressed in this form, the meaning of rj(R) becomes clear: it is the magnitude of the disk response to a unit 
perturbation in the potential. 

2.6. Self- consistent solution 
For clarity, we repeat here the important equations: 

$pcrt OR, 9) = $ (R) / port (R) cos 29 (30) 
= *%£&R,0)+*#*(R,0), (31) 

&g°(R, 0) = %&{R) %(R) cos 29, (32) 

(see eqs. Q]and[3]), and 

$*&(R, 9) = -f pcr t(R) V(R) cos 29, (33) 

(repeated from eq. I291) . The physical interpretation of these equations is that the disk response is proportional to 
the net perturbation / p0 rt, which is itself the sum of the disk response and the imposed halo perturbation. The disk 
response is opposite in sign to the halo perturbation, so /pert in the self-consistent solution must be reduced with 
respect to the imposed halo perturbation. 
The self-consistent solution can be obtained by collection equations (|3"0|) through (f3"3"| : 

MR) = /pert OR) r+^R)7^oCR) ' (34) 

In other words, the response of the disk causes the overall potential perturbation to be reduced by a factor of 1 + 
rj(R) I &o(R). All terms on the right hand side have been tabulated at the grid points Ri, resulting in a trivial evaluation 

Of /pert (R)- 

Armed with this new estimate of / pe rt, we can reexamine equation (|28p . substitute fp^ y = /pert, and calculate a 
new value of rj(R) and therefore of / per t(R)- We repeat this procedure until the maximum change between iterations 
in the quantity /pert//pert> which is a robust indicator of the relative error in / per t, is less than 10 -3 at all radii; this 
is typically achieved in 20-30 iterations. We have confirmed for some specific cases that our solution agrees to within 
~ 2% of the true solution (assumed to have converged after a very large number of iterations) at all radii, and to much 
higher precision at most radii (see Figure [1]). Adopting a stricter convergence criterion has no effect on our results. 
The relatively large number of iterations is required in order to accurately capture the sharp feature where / pC rt - > 
that is seen in the solutions (see § [3|) . 

Given / pC rt, the disk ellipticity edisk can be calculated as a function of radius directly from equation (|2"2"|) and the 
forms of the closed orbits can be calculated from equations (fTTj) - (fT4| . This provides a complete description of the 
disk. 



3. RESULTS 

In this section, we demonstrate how the disk dilutes the elliptical potential of the halo, and give examples of the net 
ellipticity induced in the disk. In § 13. 11 we demonstrate the main features of the disk-halo systems using a halo with 
constant axis ratios, while in § 13.21 we investigate how these results are affected by varying disk and halo parameters 
such as the halo concentration, axis ratio, disk scale length, and run of halo axis ratio with radius. 

3.1. Response for various disk masses 

We demonstrate the main features of our models using a fiducial triaxial NFW halo with mass M200 = 10 12 M , 
axis ratios b/a — 0.8 and c/a = 0.6, and a concentration C20 = 12. 7 These values are typical for galaxy-sized dark 
matter halos in cosmological simulations ( A llgood et al.|[2006T ). The disk rotation axis is aligned with the minor axis of 



7 M200 and C200 refer to the mass and concentration relative to the radius T200, defined such that the mean density within r200 is 200 
times the critical density. We assume the Hubble parameter h = 0.7. 



Bailin et al. 



O <D 



1 .0000 



0.1000 



£ 0.0100 



X 

o 



E 0.0010 



0.000 








20 



40 



60 



50 



100 



Iteration 



Fig. 1. — Maximum change in the solution for /pert //pert P er iteration (dashed line) and maximum difference between the solution at a 
given iteration and the true solution, assumed to have converged after 100 iterations (solid line), for the fiducial halo and 3 X 10 9 Mq disk 
of § 13.11 The vertical dotted line indicates the point where our convergence criterion is achieved. 

the halo, in agreement with the orientation of the halo angular momentum in simulations ( Bailin fc Steinmetzll2005l) . 
All disks have radial scalclcngths Rd = 3.0 kpc, with masses that range from zero up to 10 11 M©. 

Figure [2] demonstrates how the axisymmetric component of the potential, &q(R), and the rotation curve, V C (R), 
vary as a function of disk mass. The potential is given in units of V% m = G M200 /V200 • For disk masses less then 
3 x 10 9 Mq, the halo dominates the axisymmetric component of the potential at all radii, with the disk becoming 
increasingly more important with increasing disk mass beyond this. The non-axisymmetric component of the potential 
is demonstrated in Figure [3] for the 3 x 10 9 M© disk. Although the cllipticity of isopotential surfaces, e$, rises to small 
radii, the magnitude of the potential perturbation, fp^° , must vanish at small radii because Q X (R = 0) = & y (R = 0) 
while $o(i? = 0) reaches a finite value. 

Figured^ demonstrates both the initial halo perturbation (/p|J° , solid lines) and the net perturbation after including 
the self-consistent response of disk (f peI t, dashed lines). Note that even though the halo is identical in each case, /pert 
is lower for higher mass disks because the disk contributes to $ - 

IHN06I proposed that a suitable perturbation 



/port = fiso(R) ~ axe 



-x/b 



(35) 



(with x = R/r s , a — 0.1, and b — 0.098) would cause the rotation curve of a perturbed NFW profile to mimic that of a 
cored isothermal profile. We denote this as a dotted line in Figure 2^. The form of this perturbation is very different 
from the form of the perturbation that we find to be induced by a triaxial halo of uniform axis ratio, particularly once 
the self-consistent response of the disk is taken into account. 

Figure [Da demonstrates the degree to which the initial halo perturbation fp^° is diluted by the self-consistent 
response of the disk. This is equal to / P ert//pert an( i is determined from 1/(1 + rj(R)/$o(R)). The ellipticity in the 
potential vanishes in the central region where rj(R) ^> <!>o(i?). Even for a negligible disk mass, the potential in the 
innermost region is circularized. This region is larger for more realistic disks, which have a significant impact on the 
potential out to several disk scale lengths. As the disk mass increases, the form of the disk dominates both rj(R) and 
$o(i?), and therefore this function approaches an asymptotic form. 

Comparison betw een Figure HJd and the equivalent figs. 2 and 3 of IJ2KI reveals dramatically different behavior at 
small radii: in J2Kj, the "reduction factor" reaches a minimum at 1.42 Rd (= 4.26 kpc for Rd — 3.0 kpc) and then rises 
to unity, while in Figure HJd it falls monotonically to vanish at small radii. This is a direct result of the radial variation 
of /port in a physically realistic elliptical halo. Because r)(R) depends inversely on / por t (see eq. |2"51) . which must 
vanish at small radii, ?y(i?) must dominate over $q(R) in the inner regions and the potential must become completely 
circularized. 

In FigureSJ; we plot the ellipticity of the disk isodensity contours (edisk, soli d lines) and of the orbits within the disk 
(e rbit = /port 0.12, dashed lines). For a massless disk, we recover the results of IHN06I that the ellipticity rises towards 
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Fig. 2. — (a) Axisymmetric component of the potential, <t>o(^?)j for massless disks (black), and disks of mass 3 X 10 s Mq (red), 3 X 10 9 Mq 
(blue), 3 X 10 10 Mq (green), and 10 11 Mq (purple) within the fiducial halo of £j 13.11 (b) Circular velocity curve within the unperturbed 
potential. The disk contribution to the rotation curve is denoted with dashed lines. Colors are the same as in (a). 

the center of the halo. However, the presence of a massive disk changes the situation dramatically. Because /pert //port 
vanishes at small radii (Figured), the equilibrium disk is axisymmetric at small radii. Even very low mass disks, 
which contribute negligibly to $o, still respond strongly enough to the elliptical potential to cause an important change 
in the behavior at small radii. We also note that the ellipticity of disk isophotes are always greater than ellipticity of 
orbits within the disk. 

3.2. Varying disk and halo parameters 

HNSj found that the radial variation of the shape of the potential of cosmological TV-b ody halos is not consistent 
with self-similar isodensity contours (see also lJing &: Sutol 120021 ; iBailin fc Steinmet3 120051 ) . They found instead that 
the isopotential axis ratios are well fit by the following function: 



log(— or — ) = a 
a a 



tanh ( 7 log — 

r n 



(36) 



We have recomputed the self-consiste nt resp onse of disks of varying mass in a potential of this form, with the values of 
the parameters taken from halo G4 of lHNSl . whic h has a very similar mass and concentration to the halo used in § 13.11 
The results are shown in Figure [5] As noted by iHNSl , the perturbation in this case contains a peak at intermediate 
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Fig. 3. — (a) The solid curves indicate the potential along the x (lower) and y (upper) axis for the disk of mass 3 X 10 9 Mq within the 
fiducial halo of S I3.1I The dotted line indicates the axisymmetric component of the potential. The horizontal and vertical lines demonstrate 
how the ellipticity of the isopotential surfaces and the magnitude of the perturbation respectively are calculated, (b) The curves indicate 
the magnitude of the perturbation, /pert (solid), and ellipticity of the potential, e$ (dashed) for the same halo as in (a). 



radius and is much more similar to the form of /j SO required by IHN06I to fit LSB rotation curves, although unlike 
/i so the perturbation remains more prominent to large radius. However, in many cases LSB rotation curves are only 
measured out to radii of a few kpc, so it is less clear what the required form of /; so is at larger radii. The disks have 
less effect on the perturbation than in § 13.11 in particular, the radius inside which they circularize the potential is 
reduced, resulting in significantly more elliptical orbits at 1-3 kpc than for the equivalent disks in a halo with constant 
axis ratios. Disk masses higher than ~3x 10 9 M reduce the prominence of the peak in the perturbation and shift 
it to larger radii. The peak in could be moved to smaller radius, and therefore brought further into agreement 

with /iso, by reducing the r a parameter in equation (j36[) : however, this is unlikely to be a common situation, as G4 
already has by far the lowest r a value of any of the halos studied by iHNSl . 

The effect of the baryonic disk on the shape of the halo is not yet well understood. Simulations suggest that halos 
containing baryonic disks arc less elliptical th an halos composed pure ly of dark matter, and that the circularization of 
the halo occurs most strongly at th e center (Kazantz idis et ah 2004) . If the intrinsic shape of the pure dark matter 
halo is well described by the IHNSl form, which is most elliptical at the center, baryonic processes may result in a 
situation more similar to the constant axis ratio case of § 13.11 We therefore expect that the regions in which these 
models differ most strongly, 1-3 kpc, are also the regions where the unknown effect of disk formation introduces the 
most uncertainty into our models. 




R [kpc] 

Fig. 4. — (a) Magnitude of the elliptical perturbation in the potential due to the triaxiality of the halo (/JJ, solid lines), and the net 
perturbation after including the self-consistent response of the disk (/ De rt , dash ed lines). Different disk masses are indicated by different 
colors as in Figure 1. The dotted line shows the form of / pcr t proposed by HN06 to produce a rotation curve mimicking a cored isothermal 
density profile, (b) Ratio by which the input perturbation f^^° becomes diluted due to the self-consistent response of the disk, (c) 
Ellipticity of the disk isodensity contours (|ediskl> solid lines; note that edisk is negative for positive / per t in the sign convention we have 
chosen) and of orbits within the disk (e or t,it, dashed lines). 

In order to investigate how other properties of the halo and disk affect our results, we have recalculated the results 
of § 13.11 (where the halo axis ratio was assumed to be constant with radius) for the 3 x 10 9 M Q disk while varying the 
halo concentration, the b/a axis ratio, the halo mass, and the disk scale length. In more concentrated halos (FigureE]), 
the strength of the perturbation due to the halo, /port , is larger. The disk is also less able to dilute the perturbation in 
more concentrated halos. The axis ratio of the halo (Figure[7]) has a strong effect on the magnitude of the perturbation 
but has virtually no effect on the degree to which the disk dilutes the perturbation. In Figure [H] we compare halos 
of different virial mass, M2oo- In order to facilitate comparison between systems of different mass, we have kept r s 
constant by varying C200 in proportion to r2oo and kept the relative mass of the disk and halo constant. We find 
that the mass of the halo has little effect on the relative strength of the perturbation (either before or after the disk 
is taken into account), but that the resulting disk ellipticities are higher in lower mass systems. Finally, although 
the equilibrium shape of the potential is similar regardless of the disk scale length (Figure [9]) , a greater ellipticity is 
required to achieve this reduction in the perturbation for less concentrated disks, i.e. those with larger scale lengths. 

4. DISCUSSION 



4.1. Impact of Triaxial Halos on the Cusp/Core Problem 



10 



Bailin et al. 




2 



8 10 



4 6 
R [kpc] 

Fig. 5. — As in Figure[4]fbr disks in a halo potential with the same radial variation of its axis ratios as halo G4 of HNS. 

In Figure [TUl we plot the azimuthal velocities along the major and minor axes of the halo within the disk plane for 
disks of different mass in the fiducial sample halo studied in § [3l Full analysis of the velocity fields of these disks will 
be presented in a future paper. However, we note that none of these rotation curves shows the linear rise characteristic 
of a constant density co re, as e xpected given the dramatic difference between the form of / pcr t in these disks and the 
form of /i S o required by IHN06I . 



£ halo 



IS 



If we use the radial variation of the shape of the halo potential proposed by IHNSI (Figure [TTj) , for which / pcrt 
much more similar to /i so , we find that for very low disk masses the azimuthal velocity along the minor axis of the 
halo is characterized by a much more gradual rise. Observationally, such a rotation curve might be interpreted as 
indicating a constant-density core in the dark matter halo. However, for disk masses above 3 x 10 9 M Q (0.3% of the 
virial mass of the halo and just 1.8% of the baryonic mass of the system) the response of the disk removes this feature 
from the center of the rotation curve. Based on these results, we conclude that simple analyses of the shapes of halos 
are therefore not sufficient to d etermi ne whether halo triaxiality can reconcile LSB rotation curves with cuspy halo 
density profiles, as suggested bv lHN06l : full analyses that take into account the disk response are required. Preliminary 
tests on simulated velocity fields constructed usin g the r esults of this paper do suggest that triaxiality can produce 
apparent constant-density cores, in agreement with lHN06l , but a more detailed analysis includi ng many halos and lines 
of sight is needed before comparing to the observational distribution of density profile slopes (|Simon et al.ll2005l ). 

We also plot the maximum radial velocity (amplitude of non-circular motions) at each radius in Figures [10] and 111! 
Although the radial velocities are much smaller than the azimuthal velocities in all but the lowest mass disks, they are 
at a level that can be detected in observations of two dimensional velocity fields. The magnitude of the radial motions, 
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Fig. 6. — (aj Magnitude of the net elliptical perturbation in the potential (dashed lines) and the perturbation due to only the halo (solid 
lines) for disks of mass 3 X 10 9 Mq in halos with concentrations, C200 = 8 (red), 12 (blue), and 17 (green). The blue lines in these plots 
are identical to the blue lines in Figure [4] (b) Ratio by which the input perturbation becomes diluted due to the response of the disk, (c) 
Ellipticity of the disk isodensity contours (solid lines) and of orbits within the disk (dashed lines). 

which reach 5-35 km s -1 depending o n the disk mass and halo properties, are consistent with the magnitude of radial 
motions found by Si mon et al.l (|2005f ). and reproduce the observed trend for the radial motions to be negligible at 
small radii and to only become important at larger radii. 

4.2. Comparison to Observed Disk Ellipticities 

It is interesting to com pare the ellipticity of our model disks to observed values. Using two-dimensional velocity 
fields, ISimon et alj (|2005| ) found lower limits on the orbital ellipticities ranging from up to 0.175, similar to the 
orbital ellipticities in our model disks. We note that the orbital ellipticity throughout most of the disk is determined 
by the mass of the disk (Figure 0J:) , the ellipticity of the halo (Figure Efc) , and the variation of halo ellipticity with 
radius (compare Figures 0k and[5J;), with very little dependence on the concentration of either the halo (Figure [5}:) or 
the disk (Figure [9]:) , or on the global mass of the system (Figure [Sfc). We therefore predict that galaxies with large 
observed ellipticities such as NGC 4605 either lie in halos that are more triaxial than average near their center, or 
contain an unusually low fraction of their mass in their disk. The trends in intrinsic disk ellipticity tha t we p redict 
may be t e sted w ith further analysis of larger kinematic samples, such as those presented bv lGanda et alj (|2006[ ). 

iRvdenl (120061) fou nd that the distribution of isophotal shapes of galaxies in the 2MASS Large Galaxy Atlas 
(|Jarrett et al.ll2003l ), as measured in the near-infrared K s band (which is a good tracer of the stellar disk mass), 
is well fit if the intrinsic disk ellipticity distribution is a truncated Gaussian centered at 0.01 with a width of cr = 0.26. 
This corresponds to a median ellipticity of 0.18, with 68% of disks having ellipticities 0.05 < edisk < 0.37. For the disk 
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Fig. 7. — As in Figure[6]for halos with axis ratios b/a = 0.7 (red), 0.8 (blue), and 0.9 (green). 

parameters we have studied, the 18.8 mag arcsec -2 isophote at which her shapes were measured corresponds to radii 
of between 1.5 and 8.5 kpc. Our models naturally produce disk ellipticities in this range at these radii, although the 
highest ellipticities can only be produced by our least massive disks in our most elliptical potentials. 

Finally, we note that our models only take into account ellipticity in the disk induced by the dark matter halo. 
Central region s of the disk, which have high surface density and sit in an axisymmetric potential, may be unstable to 
bar formation (Berentz en fc Shlosma n 2006). This can induce additional ellipticity to the kinematic and photometric 
properties of disk galaxies (e.g.. IValenzuela et al.ll2007l ). 

5. CONCLUSIONS 

We have presented a computationally efficient method to self-c onsistently determine the dynamics of massive disks 
in triaxial dark matter halos. Our work extends the study of lJ2Kl by allowing the perturbation to the potential to vary 
with radius in an appropriate manner and by allowing the ellipticity of the disk to vary with radius self-c onsistently. 
These improvements result in qualitatively different behavior for the ellipticity of disks at small radii: IJ2KI found that 
disks counteract the halo ellipticity most strongly at 1.42 Rd and have a negligible effect at small radii; in contrast, 
we find that the effect of the disk increases monotonically to small radii, completely circularizing the potential in the 
innermost regions. 

This self-consistent radially-varying response of the disk to the halo perturbation must be taken into account when 
comparing the observed kinematic and photometric properties of galactic disks to those expected in triaxial dark 
matter halos, particularly for comparisons at small radii. When this response is calculated for plausible halo values, 
model disks have ellipticities consistent with those determined from observations of velocity fields and from isophotal 
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FlG. 8. — As in Figure l6l for halos with virial masses 

M200 = 3 X 10 11 Mq (red), 10 12 M© (blue), and 3 X 10 12 Mq (green). In order to 
facilitate the comparison, the halo scale radius r s is kept constant by varying C200 in proportion to r200i an d the ratio between the disk 
and halo mass is kept constant. 

axis ratio distributions. We also find that the radial variation of the halo axis ratios has a significant impact on the 
disk structure. Halos with axis ratios that vary with radius as suggested by cosmological simulations produce much 
more elliptical orbits in the inner disk than do halos with constant axis ratios, resulting in potential perturbations 
similar to the perturbation required to create apparent cores in galaxy density profiles. Further analysis exploring in 
detail the conditions under which core-like rotation curves might be obtained will be necessary to determine if halo 
triaxiality can resolve the cusp/core problem. 
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APPENDIX 

SECOND-ORDER TERMS AND M=4 DISTORTIONS 

We have assumed that the halo perturbation and the induced distortion in the disk are completely described by the 
m = 2 mode. This a direct consequence of only including terms linear in the small quantities / pcr t, £disk, and their 
derivatives. 
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Fig. 10. — Azimuthal velocity along the halo major axis (top lin e) an d minor axis (middle line), and the maximum radial velocity at 
each radius (bottom line) for disks of varying mass in the halo of § 13.11 NFW and isothermal rotation curves are shown for reference in 
the top right panel. 
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Fig. 11. — As in Figure [10] but for the shape of the halo potential found by HNS for their halo G4. 
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Fig. 12. — Magnitude of terms contributing to the m = 4 dist ortio n in the disk compared to £2, the magnitude of the calculated m = 2 
distortion for a disk of mass 3 X 10 9 Mq in the fiducial halo of § 13.11 

The validity of this assumption can be tested by evaluating the magnitude of the second-order terms that contribute 
to the m = 4 distortion in the disk. If we expand the potential as 



$(ft,0) = $ (R) (1 + /a (ft) cos26> + U(R) cos 40) , 



the disk surface density as 



E(JZ,( 



Sq exp 



R ( £2 OK) „„ £4 (-R) 
1- ^-^cos26»- ^-^cos46» 

Rd \ 2 2 



(Al) 
(A2) 



and include all second-order terms, then the m = 4 distortion in the disk, 64, depends on terms of order f±, / 2 62, and 
ft/2de2/dft. 8 Figure [T2l compares the magnitude of these terms to the m = 2 ellipticity for the 3 x 10 9 Mq disk in 



the fiducial halo of § 13.11 The higher-order terms are more than two orders of magnitude smaller than the first-order 
terms over most of the disk, and are also negligible in the central region where the first order terms vanish, validating 
our use of linear perturbation theory. 



As in the case of the m = 2 mode, there is also a negligible term of order Rd(fn ai^/dR. 



